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First we give a review of previous theoretical and simulation studies of 
the electrochemical double layer and then discuss our recent simulations 
and the anomalous results that we have found at low temperatures, which 
are seen in experiment but are not predicted by any existing theories. Fi- 
nally, we outline our future plans. 


1. INTRODUCTION; GOUY-CHAPMAN THEORY 


The study of the electrochemical double layer (DL) is of interest in aiding in 
our understanding of electrochemical processes, such as corrosion, at surfaces. 

Most experimental results are interpreted in terms of the Gouy-Chapman 
(GC) theory[1,2]. The GC theory is based on what is called the primitive model 
(PM). In the PM the molecular structure of the solvent, usually water, is ig- 
nored. The solvent is represented by a dielectric medium of dielectric constant 
£. In the PM the ions are represented by charged hard spheres. A further ap- 
proximation is the réstricted primitive model (RPM), where the ions are all con- 
sidered to have equal diameter,a. The GC theory goes one step further and 
considers the ions to be point charges. The neglect of the molecular nature of 
the solvent is the more severe approximation. Considering the ions to be point 
ions is not so bad and a correction called the Stern layer can be made. We will 
discuss the Stern layer shortly. 

The GC theory results from the combination of Poisson's equation, 


V’o@= -E ais) (1) 
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with the Boltzmann approximation, 


, < 0/2 


0 
a= lz qin, exp papa] , x> 02, a 


where g(x) is the electrostatic potential and q(x) is the total charge profile, q, 
is the charge and n, is the density (number of ions/volume) of an ion of species 
i,and £=1/kT , where k is the Boltzmann constant and T is the temperature. 
The parameter x is the distance from the electrode. Obviously, point ions can 
go right up to the surface, which we assume to be hard. Since the wall is 
charged, this would lead to difficulties. However, we can prevent this problem 
by not allowing the ions to get closer than o/2 to the wall. What we are doing 
is taking the ion size into account in the wall-ion interactions but neglecting it 
for the ion-ion interactions. 

The layer of thickness o/2 that is devoid of ions is called the Stern layer (or 
Helmholtz layer, or inner layer, or compact layer). It need not be regarded just 
as the distance between the ion center of a contact ion and the surface. The 
Stern layer can be considered as a solvation layer of solvent molecules at the 
electrode. There is no reason why the thickness of the Stern layer should be 0/2 
or that the dielectric in the Stern layer should have the same value as in the 
bulk. Thus, we can call d the thickness of the Stern layer and £* the value of 
the dielectric constant in the Stern layer. In principle, d and * can vary with 
electrode charge. This is often done but it is a little too empirical for this article. 

We must have electric neutrality in the bulk, 


Xa ; =0 (3) 


Furthermore, the charge per unit area in the DL must be equal in magnitude, but 
opposite in sign, to that on the electrode. Thus 
E 
Yan, 8 (x)dx =-—— =-Q (4) 
r d 4r 


The quantities Q and E are the electrode charge and the field at the electrode. 
They are related by Gauss' law. The quantity g,(x) is the distribution function 
of ions of species i. The density and charge profile of the ions of species i are 

n(x) =n,q,(x) (5) 
and 


q(x) = n.4;8 (xX) (6) 
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In the GC theory, 
g(x) =expl- fq,9(x)] (7) 


The Boltzmann relation is not valid in other theories and in simulations. 
The potential profile is given by 


p6) =-= 5 na, fit—x)e eh (8) 
i 


It is easy to verify that g(x) satisfies Poisson's equation. From Eq. 8, it is seen 
that the potential across the DL is 


rehe Be 
V=0(0)= aG ii (9) 


Equation 4 was used to obtain Eq. 9. Sometimes in simulations or in a theory, 
the ability to obtain 


V -9(0/2 j= (10) 


is used as a test of the accuracy of theory. Any theory or simulation that satis- 
fies microscopic charge neutrality, Eq. 4, will satisfy Eq. 10 automatically. 
A few moments of examination of Eq. 4 shows this. Thus, Eq. 10 is not a real 
test, at least for a simulation where charge neutrality is built in from the start. 
The right hand side of Eq. 9 may be called the Stern layer potential. The poten- 
tial g(o/2) is often called the diffuse layer potential. If a variable thickness and 
dielectric constant is used for the Stern layer, the Stern layer potential would be 
Ed/e’. 

The GC equations can be solved analytically to yield g(x), the g,(x), and V. 
For example, for a symmetric salt in the GC theory 

PAIN 


y ËL yg te? (11) 
2e fq 2 


where b= fqE/ex, x? =4nrfpq° p/e and p=} is the total density. In ob- 


taining Eq. 11 we have assumed that the ions have equal charge. Thus, q = la;|. 


It is of interest to linearize Eq. 11. Thus, at small electrode charge the GC the- 
ory gives 


y-Æ4E (12) 
2e ek 
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2. SOME RESULTS FOR THE PRIMITIVE MODEL 


In Figures 1 and 2, g(o/2) and g,(x), obtained from the GC theory, are 
compared with the computer simulations of Torrie and Valleau [3] for a mono- 
valent salt. The GC theory is based on point ions whereas the simulations are 
for ions with o = 4.25 A. The agreement is fairly good and certainly better than 
what one might expect, although it deteriorates somewhat as the electrode 
charge increases. This indicates that for monovalent salts, the size of the ions is 
not a major effect. Of course, both the GC theory and the simulations neglect 
the molecular nature of the solvent so the comparison in Figures 1 and 2 says 
nothing about the effect of the solvent molecules. 

The GC theory works less well for higher valence salts. 

0.2 
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Charge density [C/m’} 


Figure 1. Diffuse layer potential, ¥//2), for a monovalent aqueous electrolyte at 0.01, 
0.1, and 1 molar concentrations. The circles give the simulation results of Torrie and 
Valleau [3]. The dashed and solid curves give the GC and HNC results, respectively 


In Figures 3 and 4, the GC theory is compared with the simulations of Torrie 
and Valleau [4] for 1:2 aqueous electrolytes. When the monovalent ions are the 
ions attracted to the electrode (2:1 case), the results are much like the results in 
Figures 1 and 2. However, when the divalent ions (1:2 case) are the ions at- 
tracted to the electrode, the agreement is much poorer. 
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x [A] 
Figure 2. Normalized density profiles, g (x)=n,(x)/n;(%), for a 1 molar monovalent 


aqueous electrolyte. The circles and curves have the same meaning as in Figure 1 
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Figure 3. Diffuse layer potential, y(o/2), for an asymmetric aqueous electrolyte at 


0.005, 0.05, and 0.496 molar concentrations. The circles give the simulation results of 
Torrie and Valleau [4]. The curves have the same meaning as in Figure 1 
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x [A] 
Figure 4. Normalized density profiles, g(x) =n,(x)m,(), for a 0.496 molar asymmetric 
electrolyte. The circles and the curves have the same meaning as in Figure 3 


In addition to the GC theory, the hypernetted chain (HNC) approximation 
and the mean spherical approximation (MSA) have been applied to study DLs 
in the primitive model. The MSA can be regarded as a linearized version of the 
HNC approximation. The advantage of both is that the size of the ions can be 
considered in a self consistent way. They have the further advantage that they 
can be applied to nonprimitive models, where the molecular models of the sol- 
vent can be considered. 

Blum [5] has applied the MSA to study DLs in the PM. He obtained 

E 


= ——_. 13 
V OT) (13) 


where J" is a renormalized screening parameter that is related to x by 
k =2I(1+TIo0) (14) 


Expansion of Eq. 13 in powers of x leads to Eq. 12 as one would hope. 
As is seen in Figures 1-4, the HNC results are in better agreement with the 
simulation results than are the GC results. The HNC results fail to agree with 
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the simulation results at high electrode charge. Density functional (DF) theory 
[6] and the HNC2 or MSA2 theory [7] give excellent agreement. 


3. MOLTEN SALTS; LOW TEMPERATURE ANOMALIES 


One promising application of the PM is the study of DLs in molten salts. The 
advantage of molten salts is that there is no solvent so the PM should be a good 
model. The disadvantage of molten salts is that, because the dielectric constant 
is 1 or near 1, the value of the coupling constant, fe*/o , is very large. This 
makes simulations and theory very difficult. 

One puzzling feature of molten salt DLs is that according to experiments the 
capacitance 


Q 

C s (15) 
in such systems increases with increasing temperature. The GC theory predicts 
that 

dinC 1 

d\nT 2 (16) 
which has the opposite sign. Furthermore, every theory that has been applied to 
this problem yields a negative temperature derivative. Thus, one is faced with 
the question is to whether the seemingly anomalous temperature dependence of 
the experimental capacitance is due to any failings of the PM or to the deficien- 
cies of the theories. 

To answer this question, we have simulated the molten salt DL using the 
PM. As is seen in Figure 5, the slope of C* vs. T* is negative at high tempera- 
tures in qualitative agreement with the GC theory but becomes positive at low 
temperatures (the definitions of reduced quantities are found in the captions of 
the figures). Since molten salts exist at high temperature, it may seem odd to 
compare the molten salt behaviour with that of the simulations at low T*. How- 
ever, because ¢ =1 in molten salts, they are low T* systems. 

Before this work began we thought that the temperature behaviour of molten 
salt DLs was due to the fact that the density was high. From Figure 5 it is ap- 
parent that the anomalies are because of the low temperature of the system. This 
suggest that the temperature dependence of the capacitance of aqueous electro- 
lytes might exhibit the same positive slope at temperatures below room tem- 
perature. As is seen in Figure 6, at low temperatures DLs in aqueous electro- 
lytes exhibit the same behaviour as do molten salt DLs. Since starting this 
study, W. Schmickler has drawn our attention to the study of Hamelin et al. [8] 
who found that the DL capacitance of HClO, - 5.5H2O at temperatures between 
200 K and 300 K decreases with decreasing temperature. 
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Figure 5. Reduced capacitance, C* =Qa/V, for a model salt double layer as a function 
of the reduced temperature, T° =okT/e?. The reduced density is typical of a liquid: 
p` =po’ =0.64. The dotted line is to guide the eye 
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Figure 6. Reduced capacitance, C* = Qo/V , for the double layer in the PM as a function 
of the reduced temperature, T° = ok7/e*, near the critical temperature. The reduced den- 
sity is low: p° = po? =0.04 
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Regrettably, the MSA cannot account for the positive slope of C* vs. T° 
that we have seen in our simulations. The HNC approximation has not been 
tried in this application but in view of its relation to the MSA there is no reason 
to expect it to be successful. Disappointingly, the DF theory fails also. The 
HNC2 and MSA2 approaches have not been used yet for this problem. 
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0 0.02 0.04 0.06 0.08 
p* 
Figure 7. The adsorption isotherms of the PM as a function of the reduced density at 
two temperatures close to the critical temperature. The electrode charge density is 
Qc’ /e = 0.03. The dotted lines are to guide the eye 


According to simulations, there is a critical point in electrolytes at T* ~ 0.05 
and p` ~ 0.025 [9, 10, 11]. The value of 7” is better established than the value 
of p`. In any case, Henderson et al. [12] have presented qualitative arguments 
that the adsorption isotherm, 


r, = Yn Jle(a)-thx (17) 


would be large in the near vicinity of the critical point. As is seen in Figure 7, 
this is supported by our simulation studies. 


4. NONPRIMITIVE MODELS 
As we have mentioned, the simplest method for taking the solvent into ac- 


count is to give the Stern layer a different thickness or more importantly a low 
dielectric constant. This latter semiempirical change mimics the presumed fact 
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that the solvent molecules are oriented near the electrode. Using the latter modi- 
fication, Eq. 11 would be replaced by 


EG, À b 
V =—~-+—sinh" — 18 
2e Bq j) ajap 
and Eq. 12 would be replaced by 
EA (19) 
2e & 


The main consequence of Eqs. 18 and 19 is that the reciprocal of capacitance 
would be the sum of two terms 


eE T ET (20) 


The first term is the Stern layer term and the second term is called the diffuse 
layer term and is given by the GC theory. If Eq. 20 makes any sense, we can 
plot C™' versus C;) and get a straight line of slope 1. As is seen Figure 8, the 
experimental capacitances [13] do fall on straight lines to a good approxima- 
tion. The values of C™ that are plotted in Figure 8 are differential capacitances, 
dQ/dV , rather than the integral capacitances, Q/V . The distinction is not mate- 
rial to our discussion. 

The value of C™ at high concentrations (1/C,,, =0) give the values of C;} 
as a function of electrode charge. These values are plotted in Figure 9. Clearly, 
if Eq. 18 were to be used, ¢° would have to be a function of electrode charge. 
We will not do this. Our attention will be directed to the zero charge case. 

The next step in sophistication would be to represent the solvent by hard 
spheres and to retain the dielectric background. The hard sphere solvent would 
account for the excluded volume effect. This has been done by Henderson and 
Lozada-Cassou [14] and by Tang et al. [15]. The latter authors called this the 
solvent primitive (SP) model. Tang et al. have made simulations for this model. 
At zero electrode charge, the MSA gives Eq. 12 or 19. 

The main result for the SP model is to cause strong oscillations in the density 
profiles. Since they affect the positive and negative ions in much the same man- 
ner, the potential and capacitance are not so strongly affected. Even so if the 
profiles differ appreciably from the GC predictions, the microscopic picture is 
quite different. Oscillations in profiles have been seen in experimental determi- 
nations of the force between macroscopic bodies immersed in electrolytes [16]. 
Since electrochemical reactions are, presumably, affected by the microscopic 
situations, many ideas in electrochemistry should be reexamined. 
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Figure 8. The inverse of the experimental differential capacitance as a function of the 
inverse of the diffuse layer capacitance, as given by the GC theory. The symbols are the 
experimental results of Parsons and Zobel [13] for an aqueous solution of NaH,PO, at 
25 °C near a hanging mercury electrode. The dotted lines give straight lines of slope 1. 
The solid line gives the full MSA result for ion-dipole electrolyte 
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Figure 9. Plot of the intercepts of the straight lines in Figure 8 with the C7} =0 axis as a 
function of the electrode charge density 
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We can go one step further and give the hard sphere solvent molecules a di- 
pole moment, yw, and, thereby, eliminate the need for a dielectric background. 
Carnie and Chan [17] and Blum and Henderson [18] have solved the MSA 
equations for this ion-dipole model. As one would expect, the profiles are os- 
cillatory. Instead of Eq. 19, they obtained an explicit but complex result, which 
upon expansion in powers of x , gives 


Vas ey SS... (21) 
ex 2e A 
where À is related to € by 
V7(1+A)* =16e (22) 


For ¢=1 (molten salts), 4=1 and for ¢=78 (aqueous solutions), 4 = 2.65. 
Thus, À is a weak function of € . 
Comparison with Eq. 19 yields 


3 are (23) 


which gives a value for £* that is similar to that used in fits to experiment. 


sim. (our results) 
* sim. (Berendsen et al.) 
experimental 
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r [nm] 
Figure 10. The radial distribution function of water obtained from simulation using the 
SPC/E model, compared to the simulation results of Berendsen et al. [20] and the ex- 
perimental curve 
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This does not mean that the MSA validates the idea of an inner or Stern layer 
in series with a diffuse layer, even though Eq. 21 looks similar to what would 
result from such a picture. The potential results from solvent and ion contribu- 
tions, each of which comes from an integration over all space. Thus, the solvent 
and ionic density profiles extend over all space. There is no inner layer where 
the solvent dominates. Furthermore, £* is merely a parameter and is not a di- 
electric constant in any part of the space. It is our view that the dielectric con- 
stant is a microscopic quantity and it is not meaningful to talk about the dielec- 
tric constant in a macroscopic region such as the Stern layer. 

We have plotted the full MSA inverse capacitance in Figure 8. At low con- 


centrations (larger C}, ), the MSA curve is a straight line of slope 1. However, 


at high concentrations, the MSA inverse capacitance falls below the straight 
line. Interestingly, so do the experimental values. 

We have made some simulations [19] of the ion-dipole model of the DL. We 
are developing a simulation code for the DL using a realistic model of water. As 
of the date of the writing of this article, we do not have results for the DL with 
this realistic solvent. However, a water-water correlation function for bulk wa- 
ter is plotted in Figure 9. We hope to adopt our program to the DL geometry 
and be able to report results in the near future. 
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